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y—( • Abstract. We analyze the Parkes quarter-Jansky flat-spectrum sample of QSOs in terms of space density, including 

00 ' the redshift distribution, the radio luminosity function, and the evidence for a redshift cutoff. With regard to the 

>^ , luminosity function, we note the strong evolution in space density from the present day to epochs corresponding 

' to redshifts ~ 1. We draw attention to a selection effect due to spread in spectral shape that may have misled 

I other investigators to consider the apparent similarities in shape of luminosity functions in different redshift 

. shells as evidence for luminosity evolution. To examine the evolution at redshifts beyond 3, we develop a model- 

CIh' independent method based on the Knax test using each object to predict expectation densities beyond z — 3. 

^ , With this we show that a diminution in space density at 2 > 3 is present at a significance level > 4cr. We identify 

t I ■ a severe bias in such determinations from using flux-density measurements at epochs signiflcantly later than that 

' of the finding survey. The form of the diminution is estimated, and is shown to be very similar to that found for 

, QSOs selected in X-ray and optical wavebands. The diminution is also compared with the current estimates of 

II" ' star-formation evolution, with less conclusive results. In summary we suggest that the reionization epoch is little 

I infiuenced by powerful fiat-spectrum QSOs, and that dust obscuration does not play a major role in our view of 

■ the QSO population selected at radio, optical or X-ray wavelengths. 
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1. Introduction 

This is the last in a series of three papers describing the 
results of a program to search for high-redshift radio-loud 
QSOs and to study the evolution of the flat-spectrum QSO 

population. 

Paper 1 l)jackson et al.ll2002() set out the sample, dis- 
cussing selection, identification and reconfirmation pro- 
grammes to determine the optical counte rparts to the ra- 
dio sources. Paper 2 l|Hook et alJ l2003l) presented new 
spectroscopic observations and redshift determinations. 
This paper considers the radio-loud QSO space distribu- 
tion, the epoch-dependent luminosity function, the evi- 
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dence for a redshift cutoff provided by the sample, and 
the form of this cutoff. 

Paper 1 described how the identification programme 
for 878 flat-spectrum radio sources selected from the 
Parkes catalogues yielded a near-complete set of optical 
counterparts. Indeed for the sub-sample at declinations 
above —40° with flux densities above catalogue complete- 
ness limits, only one source remains unidentified. Of the 
379 QSOs in this sub-sample, 355 have measured redshifts, 
obtained from earlier observations and the extensive spec- 
troscopy programme described in Paper 2. This relative 
completeness is ideal for studies of space density, as it be- 
comes possible to map the entire "quasar epoch" with a 
single homogeneous sample, having no optical magnitude 
limit and free of obscuration effects. In fact a sub-samp le of 
object s from an earlier analysis was used by [Shaver e t alJ 
l|l996j) to study the evolution of QSO space density at 
high redshifts. The study showed that the space density 
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of high-luminosity radio QSOs decreased significantly at 
redshifts beyond 3. Pr eliminary data were also used by 
iJackson fc Walll lll999l) in considering a dual-population 
scheme of space densities for unified models of QSOs and 
radio galaxies. 

General features of the luminosity function and its red- 
shift dependence haye long been established for QSOs 
selected at optical and radio waveleng ths (e.g. Longair 
llflfifil ISchmidtJlTflfii iFanti et a,ljfl97.-^ . Powerful evolu- 
tion is required, similar in magnitude for selection at 
either waveband; the space density of the more lumi- 
nous QSOs at redshifts of 1 to 2 is at least 10^ that 
at the present epoch. It has been hotly debated as to 
whet her the form of th is change is luminosity evolution 
(e.g. iBovle et al.lll988() or luminosity-d ependent density 
evolution fe.g.bunl op & Peac ock"l990V It does not mat- 
ter: physical models are not available that require either 
form, although it is clear that lumino sity evolution re- 
sults in lifetimes of non-physical length iHaehn elt fc Reei 
ll99,Sl and references therein) . The space density of radio- 
selected QSOs, constituting some 10 per cent of aU QSOs, 
generally appe ars to par allel th at of optically- selected 
QSOs fe.o. Schmidt et alj.1991: Stern et alJl200ol) . 

There are many reports of a redshift cutoff in the lit- 
erature: paper after paper speaks of 'the quasar epoch', 
'a strongly-evolving population peaking at a redshift of 
about 2', or 'the quasar redshift cutoff' without spe- 
cific reference. For optically-selected QSOs, several clas- 
sic studies demonstrated that such a cutoff does exist 
lISchTTiidt et a,1 .1119911 IWa,rren et a,l.iri994llKeTinefick et a,l.l 
Il995(l . The Sloan Digital Sky Survey (SDSS) has now 
found QSOs o ut to redshifts beyo nd 6, and analyses of the 
space density (jPan et al.ll2001al3 lbh provide the strongest 
evidence to date of the drop in space density beyond 
2 = 3. X-ray surveys now appear to show that the X-ray 
QSO population exhibits a decline at high reds hifts sim- 
ilar to that found for optically-se lected QSOs l)Hasingeil 
^003i Bargcr et al. 2003; Silverm an et al.ll2004|) . But do 
all these observations indicate a real diminution or - at 
least at optical wavelengths - could it be d ue to a dust 
screen (Hcislcr & Ostrikcr 1988; Fall & Peil ll993l) ? It is 
here that radio-selected samples such as the present one 
can provide a powerful check: if a significant diminution 
is seen in the radio luminosity function, it cannot be the 
result of dust obscuration. iDunlop fc Peacock! l)l99Cil pre- 
sented some evidence for just such a cutoff of the ra- 
dio luminosity function (RLF) for fiat-spectrum (QSO- 
dominated) populations; and a n earlier analysis of a sub- 
sample from the present wo rk l|S haver et"al 1ll99()t) added 
confirmation. More recently IVigotti et''al!r|2'003(l defined 
a complete sample of 13 radio QSOs at z ^ 4, from which 
they concluded that the space density of radio QSOs is 
a factor of 1.9 ± 0. 7 smaller than tha t of si milar QSOs 
at z ~ 2. However, Ijarvis fc Rawlingi l|200d|) questioned 
these radio-QSO results, focussing on the possible effects 
of spectral curvature. 

A possible dust screen has serious implications for 
the interpretation of the Hubble diagram for SN la su- 



pernovae. Assuming no obscuration, current results from 
the SCP (Supern ova Cosmology Pro ject) collaboration 
llKnop et al.ll2003|) and the Hi-z team jTonrv et al"ll2003|) 
favour an f2m = 0.3, — 0.7 universe. Two further 
related issues make delineation of the QSO epoch very 
important: galaxy formation, and the reionization of the 
Universe. 



1.1. Galaxy formation 

The dramatic cosmic evolution of radio galaxies and QSOs 
stood as a curiosity on its own for over 30 years since the 
birth of the idea (Rylc 1 955.1. clouded as it was in the 
source-count controversy llScheueT<IT99(i . It is relatively 
recently that corresponding evolution has been deline ated 
for th e star- formation rate (Lilly et al. 1996: Madau et alJ 
I199(tI) and for galaxy evolution, particularly blue galaxies 
llEllislll999t) . The correlation between star-formation rate 
and AGN space densit y (Wall 1998) strongly s uggested 
a physical connection l|Bovle & TerlevichI il99a|) . Before 
the emergence of the Lilly-Madau plot of star-formation 
history, it was recognized that the model of hierarchi- 
cal galaxy development in a Cold Dark 
Universe would result in a 'quasar epoch' 
iHachrielt fc R,ees .1993,') . The issue of 'quasar epoch' and 
'redshift cutoff' has therefore assumed particular impor- 
tance in consideration of galaxy formation in low-density 
CDM universes. The very existence of any high-redshift 
QSOs sets constraints on the epoch of formation of the 
first galaxies. Haehnelt (1993) showed how the then-new 
COBE normalization (Smoot et al. 1992) together with 
the QSO luminosity f unctio n at high redshifts as mea- 
sured by iBovle et all lll99ll) . provided substantial infor- 
mation on the initial fluctuation spectrum and the matter 
mix. He found that the z = 4 luminosity function ex- 
cluded an initial-spectrum index of n < 0.75 or a Hot 
Dark Matter fraction > 25 per cent. Relevant to the cur- 
rent view of the low-ma tter-density CDM Unive rse, he 
found that f^A < 0.75. iHaehnelt fc Reed l|l993|) devel- 
oped a model for the evolution of the QSO population 
based on the existence of ^ 100 generations and link- 
ing the QSO phenomenon with the hierarchical build-up 
of structure in the Universe. The evolution of host ob- 
jects is mirrored in the evolution of the mass of newly 
formed black holes; only a moderate efficiency for forma- 
tion of an average black hole is necessary to model the 
luminosity function. The model suggested that nearly all 
galaxies are likely to have passed through a QSO phase. 
[Kauffmann & Haehnelt (2000) produced a more sophis- 
ticated model by incorporating a simple scheme for the 
growth of supermassive black holes into the CDM semi- 
analytic models that chart the formation and evolution 
of galaxies. In addition to reproducing the observed re- 
lation between bulge luminosity and b lack-hole mass in 
nearby galaxies (^agorrian et abl l 19981) . the model is able 
to mimic the enormous increase in the QSO population 
from redshift to 2, a feature that the Hachnelt-Rees 
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model was able to describe only qualitatively. Their con- 
clusion: "Our results strongly suggest that the evolution 
of supermassive black holes, quasars and starburst galax- 
ies is inextricably linked to the hierarchical build-up of 
galaxies." 

1.2. Reionization 

The paradigm of hierarchical structure growth in a 
CDM universe has long suggested that after the re- 
combination epoch a,t z 1500, the reionization of 
the Universe took place at r edshifts between 6 and 20 
(e.g. iGnedin fc Ostrikeilll997t) . Th is reionization is pre- 
dicte d to be patchy and gradual I Miralda-Escude et al.1 
although some i nodels indicate that it should hap- 
pen q uite rapidly (e.g. ICen fc McDonaldll2002t iFan et all 
l2002h . Two major observational advances support the 
'patchy and gradual' scenario. Firstly, SPS S discovery of 
QSOs at redshifts of 6 or more ijPan et alllionnl l2001at 
iBecker et aO 1200 1|) has given a glimpse of what may 
be t he end of the epoch of r eionization: the first com- 
plete iGunn fc PetersonI trough has b een observed 
in the z = 6.29 QSO SP SS 1030-^0524 llBecker et all 
l200l[|Pentericci et al.l20o3) and a second has been seen in 
z = 6.43 QSO SPSS J1148+5251 JFa,n et a,1.ll2003t) . There 
is disagreement as to whether this ma rks the true end 
of reionization llSongaila fc Cowiell2002^ : but the sugges- 
tion is that it may be essentially complete by z ~ 6 — 7. 
Secondly, the detection of polarized anisotropies with the 
Wilkinson Microwave Anisotropy Probe (WMAP) has re- 
sulted in a measurement of t he optica l depth r ^ 0.17 
to Th ompson scattering ( Bennett et al. >2003-: Kogut et al.l 
|2003() . implying a reionization redshift of 17±5. The CMB 
is sensitive to the onset of ionization, while Gunn-Peterson 
troughs are sensitive to the late stages, the cleanup of 
remaining HI atoms. Resolving the large uncertainties 
in these redshifts could yet result in a rapid reioniza- 
tion scenario. Neverthe less several recent papers (see e.g. 
iHaiman fc Holderl2003l) address the complex and interact- 
ing suite of physical mechanisms that may be at play dur- 
ing an extended 'patchy and gradual' reionization epoch 
6 < z < 17. 

In either a fast or a gradual scenario, identifying the 
source of this reionization as well as epoch is of vital im- 
portance for such interconnected reasons as: 

— The role of reionization in allowing protogalactic ob- 
jects to cool into stars, 

— The small-scale temperature fluctuations in the CMB 
and how these are influenced by patchiness in the 
reionization, and 

— The epoch of the first generation of stars, or galaxies, 
or collapsed black-hole systems. 

It is most likely that the reionization is via photoion- 
ization by UV radiation from stars or QSOs, rather than 
collisional ionization in e.g. blas t waves from t he explosive 
deaths of Population III stars l|Ma,daiill200(i . It may be 



possible to detect this reionization epoch directly as a step 
in the background radiation at radio frequencies between 
70 and 240 MHz (redshifted 21-cm H I) or in the infrared , 
0.7 to 2.6 fim, from H recombination (|Shaver et alJll999ll . 
It may be possible with results from the Planck mission 
to identify features in the CMB that identify what the 
predominant mechanisms are; and it may be possible to 
detect the UV sources responsible for the ionizing flux 
at z ^ 10 — 20 with th e James Webb Space Telescope 
l|Haima,n fc HoldeHEooi. 

QSOs have long been prime candidates for this 
reionization. However the apparent decline in space 
density (f rom the evidence summarized above and by 
iMadau et al. 1999), is inconsistent with this interpreta- 
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tion. IMadaia ((2000) showed that in the face of this appar- 
ent diminution, UV luminosity functions of Lyman-break 
galaxies (LBG) provide 4 times the estimated QSO con- 
tribution at z = 3. It is now commonly accepted that 
such objects (or their progenitor components) take on the 
mantle. The formation of short-lived massive stars in such 

H l200lh . 



galaxies provides the UV photons IjHaehnelt 
although QSOs may supply a significant fraction of the 
UV background at lower redshifts. 

Because the cooling time is long, the low-density IGM 
retains some memory of when and how it was ionized. 
Several investigators h ave found a peak in temperature 
oOhe IGM at z ~ 3 jSchave et a,lj 1200(1 iTheuns et a,lJ 
|2002|) close to the peak of the 'quasar epoch'. Moreover, 
observations of several QSOs at the wavelength of Hell 
Lya near z — 3 suggest delayed reioniz ation of He I, wit h 
the process not yet complete by z = 3 iKriss et al .11200 ll) . 
The implication is that the QSO ionizing photons coinci- 
dent with the peak in activity both rcionize Hel and dump 
entropy into the IGM to raise its temperature. 

In all of these aspects, it is clear that conclusions on 
ionizing flux from QSOs are dependent on poorly de- 
termined high-power regions of luminosity functions, on 
apparent cutoffs observed primarily in optically-selected 
samples, and then only for the most luminous QSOs. 

It is a primary purpose of this paper to determine the 
radio luminosity function using the near-complete data 
of the present sample, and to examine the evidence for 
a redshift cutoff. Before this, we discuss the populations 
involved in the flat-spectrum sample by examining the 
N{z) relation (§0). Subsequent to the RLF determination 
in §0 we consider the issue of a redshift cutoff (§0)), and 
the form of this cutoff. In § El we construct an overall 
picture of epoch dependence of space density for radio- 
loud QSOs. We compare this with the parallel results for 
QSOs selected at optical and X-ray wavelengths, and with 
the behaviour of star-formation rate with epoch. The final 
section (§ summarizes results from this paper and our 
preceding two papers. 

2. The redshift distribution 

For a sample of objects complete to some flux-density 
limit, the redshift distribution, N(z), gives preliminary in- 
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formation on the epoch of the objects, and allows the most 
direct comparison with other samples. The redshift distri- 
bution gives direct information on neither the luminosity 
function nor its epoch dependence; however it provides es- 
sential data for use with other data such as source counts 
to enable the construction of epoch-dependent luminosity 
functions. There have been many versions of this. Most 
are a variant on either the T4iax method (jSchmidt,.1968tl 
or the technique of defining the luminosity distribution 
ijWall et al.lll98f]l:IWal]l ll98,'^'>. obtained when a complete 
N{z) is available at one flux-density level at least. 

Such modelling processes now make use of statisti- 
cal techniques to incorporate data sets of varying com- 
pleteness at many frequencies and flux-density levels. The 
sample described here rep resents only one such data set, 
more complete than most. IPunlop fc Peacocl] ((199Q') car- 
ried out the most extensive such modelling. They took 
as a starting point two populations, 'flat- spectrum' and 
'steep-spectrum' radio sources, now broadly considered in 
the light of unified models as beamed radio sources (radio- 
loud QSOs and BLLac objects) and their unbeamed pro- 
genitors, or hosts (FRI and FRII radio galaxies). All these 
objects are deemed to be powered by accretion-disk / ro- 
tating black-hole systems from which a pair of opposing 
relativistic jets feed double radio lobes. The single axis 
is collimated during the feeding process by rotation of 
the black-hole system. The beamed objects, QSOs and 
BL Lacs, beamed because of relativistic ejections of com- 
ponents along axes aligned with the line-of-sight, have ra- 
dio structures dominated by relativistically-boosted core 
emission. This core emission shows the effects of syn- 
chrotron self-absorption and therefore has a flat or in- 
verted radio spectrum. The radio emission from powerful 
radio AGN whose axes are not aligned with the line-of- 
sight is dominated by their steep-spectrum lobe emission, 
on the large scales of lO's of Kpc up to lOO's of Mpc. 
The dichotomy between beamed and unbeamed objects 
as evidenced by their integrated radio spectra is shown in 
Fig. [2 

The widely-used Dunlop-Peacock models of the lumi- 
nosity functions may be simply tested against the present 
data by means of the N{z) distributions that they predict. 

We constructed a redshift distribution from the sam- 
ple of Paper 1 as follows. We selected all sources with 
5*2.7 GHz > 0.25 Jy in regions for which the 2.7-GHz flux- 
density limit was 0.25 Jy or less, and at declinations 
2.5° >S> -40°. We refer to this as Sample 1 and the to- 
tal area it covers (Fig. 1 of Paper 1) is 2.676 sr. The source 
composition, identification and redshift data for this sam- 
ple are shown in Table Choice of the declination limit 
comes from both identification and radio-spectral com- 
pleteness; see §3. 

The entries in the identification column. Table ^ re- 
fer to (QSO)s, (BL) Lac objects, (G)alaxies, (Obsc)ured 
fields, and (e) not identified for reasons discussed in 
Paper 1. As reasonable approximations, the 20 QSOs with- 
out measured redshifts were assumed to have the same 
redshift distribution as those with redshifts; likewise the 
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Spectral index 2.7 - 5.0 GHz 

Fig. 1. The spectral index distribution for the compila- 
tion of sources from Parkes surveys known as PKSCat90. 
The sub-sample selected has <S'2.7GHz > 0.25 Jy, in re- 
gions surveyed to a limiting fiux density of 0.25 Jy or less. 
The hatched area shows the sources identified as QSOs 
or BLLac objects. These beamed objects dominate the 
'flat-spectrum' region of the diagram. 

Table 1. Sample 1: an all-source sample selected from 
the source list of Paper 1 in order to estimate the redshift 
distribution: 52.7 ghz > 0.25 Jy, < 0.25 Jy, 2.5° > 
6 > -40°. 



Ident 


No Redshift 


Redshift 


Total 


QSO 


20 


308 


328 


BL 


34 


9 


43 


G 


57 


27 


84 


Obsc 


2 





2 


e 


1 





1 


Totals 


114 


344 


458 



unmeasured redshifts of the 34 BLLac objects were as- 
sumed to have the same distribution as those measured. 
Such an approximation is inappropriate for the galaxies, 
however. A crude Hubble diagram was plotted for the 27 
galaxies with redshifts and a simple polynomial was fltted 
to make rough estimates of the redshifts for the remain- 
ing 57 galaxies based on their B magnitudes. Finally the 
3 Obsc and e sources were assumed to have the same red- 
shift distribution as the total sample; the N{z) obtained 
by adding the QSO, BLLac and galaxy redshifts was sim- 
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ply scaled by (344 + 111 + 3)/ (344 + 111) to obtain the 
final N{z) of Fig. 




1 2 3 4 5 

Redshift 



Fig. 2. The redshift distribution (histogram) for the 
sources of Sample 1 (Table^. The hatched area shows the 
redshift distribution for beamed objects alone, the QSOs 
+ BLLac objects, while the clear region represents the 
galaxies. The 6 dotted lines show the appropriately-scaled 
distributions predicted by the flat-spectrum (a > —0.5) 
components of the Dunlop-Peacock (1990) luminosity- 
function models, with the dashed line distinguishing their 
pure-luminosity-evolution model. The solid line represents 
the average of 6 of the models, omitting the model showing 
the very steep rise to z 5. The symbol -I- line systems 
show the predictions from the dual-population models of 
Jackson & Wall (1999), red representing all beamed ob- 
jects (QSOs + BLLacs), blue for QSOs only. 



Dunlop & Peacock derived luminosity functions for 
their two-population model, flat-spectrum and steep- 
spectrum radio sources, representing the luminosity func- 
tions as polynomials over the surface (p, Pradio, z), and ob- 
taining coefficients by best-fitting to multi-frequency sur- 
vey data including source counts and redshifts. Different 
models resulted from different starting points and fac- 
torizations of the epoch-dependent luminosity function. 
Their division between flat-spectrum and steep-spectrum 
sub-populations took spectral index a = —0.5 as the 
dividing criterion. The predictions of redshift distribu- 
tions from the flat-spectrum portions of Dunlop-Peacock 
luminosity functions are shown in Fig. |21 In order to 
scale these to our spectral-selection criterion, we used the 
spectral-index histogram of Fig. ^ the ratio of objects 
with a > —0.4 (our selection criterion. Paper 1) to those 
with a > -0.5 is 1060/1275 = 0.831. 



In view of uncertainties in spectral index and of equat- 
ing the flat-spectrum population of Dunlop & Peacock 
with compact radio sources, the overall agreement is good. 
The form of the decline in N{z) to higher redshifts is im- 
pressively described by the Dunlop-Peacock models. Two 
models stand out in Fig. |21 One model with a space- 
density cutoff at z = 5 predicts a redshift distribution 
greatly at variance with observations, showing a dominant 
spike in the distribution at redshifts just below this cutoff. 
It has been left out of the averaging process. The pure- 
luminosity-evolution model, shown as the dashed line, is 
distinct in having a quicker rise and flatter maximum than 
the others. These two features provide a better represen- 
tation of the data in the range < z < 1.5 than do the 
other models. 

The good fit of the Dunlop-Peacock models to the to- 
tal N{z) distribution for flat-spectrum objects does not 
imply a good description of the N(z) for beamed ob- 
jects (hatched area. Fig. ^ alone. The Dunlop-Peacock 
models clearly rely on the presence of low-luminosity 
flat-spectrum galaxies for the quality of overall fit; the 
'flat-spectrum' models describe the beamed objects alone 
rather poorly. 

Models considering populations in terms of beamed 



and host object 


were 


developed by Jackson and Wall 


(WaU & Jackson 


,1997; 


Jackson & Wafl 1999). The N{z) 



predictions from these models are shown in Fig. 13 
Agreement is reasonable; normalization is correct, and the 
forms of the curves are similar. This agreement is expected 
on the basis of the fit of the model to the 5-GIIz source 
count and the incorporation of a redshift cutoff in the 
model evolution. The models over-predict objects at z > 2, 
due primarily to a lack of constraint on the evolution of 
low-luminosity sources. 

3. The Radio Luminosity Function (RLF) 

Completeness of identifications enables the radio luminos- 
ity function to be construct ed in a straigh t forward way , 
using the l/V„h, approach fechmidtl Il96i iFeltenI [l976t 
^vni & BahcaU 1980). The contribution of each object to 
space density is calculated as the reciprocal of the observ- 
able volume, the volume defined by the redshift range(s) 
in which the object can be seen. Because the sample is op- 
tically complete, only radio data (apart from the redshifts) 
are relevant in defining this range. 

An appropriate sub-sample for this calculation is that 
referred to as Sample 2 in Table [3 Selected from the 
catalogue of Paper 1, it includes all the QSO identifica- 
tions with flux densities above survey limits and within 
the declination range -1-2.5° to —40°. Defining Vots re- 
quires knowledge of the radio spectrum both above and 
below the survey frequency. Above 2.7 GHz, there are 
the 5.0-GHz data of the Parkes catalogues for all sources 
in the 2.7-GHz surve ys, fltuc densities a t 8.4 GHz for 
many of these sources llWright et al.lll99nh . and about 40 
8.87-GHz flux densities for some of the brighter sources 
i Shimmins fc Wal]lll973l) . Below 2.7 GHz, flux densities 
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exist for most me mbers of Sa mple 2 at 365 MHz from 
the Texas survey llDoudas et al...l99ffl. and at 1.4 GH z 
from the NRAO VLA sky survey l|Condon et alJll998() . 
The Texas survey covers the sky at declinations down 
to -35.5° and the NVSS down to -40°. As a compro- 
mise between sample size and spectral completeness, the 
sub-sample chosen for definition of the RLF, Sample 2 of 
Table El was therefore taken to have a southern declina- 
tion limit of —40°. Most of the area surveyed at 2.7 GHz 
in this range has a completeness limit of S'2.7ghz = 0.25 
Jy, but some regions have limits of 0.10, 0.20 and 0.60 Jy; 
see Fig. 1 of Paper 1. 

Table 2. QSO samples for RLF and redshift-cutoff anal- 
yses. 



Sub-sample Sample 2' Sample 3 



with measured z UKST ID 342 242 

with measured z CCD ID 13 10 

total QSOs with measured z 355 252 

no measured z UKST ID 17 11 

no measured z CCD ID 6 4 

total QSOs with no measured z 23 15 

no ID (no measured z) 1111 



Totals QSOs + non-ID 379 379 



These data are displayed in Fig. O First impressions 
are that the curves slide sideways, suggesting simple lumi- 
nosity evolution, as deduced from similar beha viour in red- 
shift shells for optical luminosity functions ijBovle et alJ 
1988). However, the transition from the lowest redshift 
range {z < 0.5) to the next redshift range (0.5 < z < 1.0) 
is not described by a lateral shift. It is possible that the 
lower bin is contaminated by unbeamed objects such as 
Seyfert galaxies and elliptical-galaxy cores; such objects 
may have entirely different central engines and different 
luminosity functions as a result. The data of Fig. |31 shown 
in integral form (Fig. ^ suggest that the RLF changes in 
form right out to the 1 < z < 2 shell, and it is improbable 
that contamination by unbeamed objects persists beyond 
z > 0.5. A closer investigation by both morphology and 
spectrum is needed to determine if removal of unbeamed 
objects could 'save' luminosity evolution. However, it is 
probably beyond saving. For example, from QSOs discov- 
ered in the SDSS survey, iFan et all l|2001( j) noted that the 
high-power end of the QSO luminosity function appears 
flatter than that at lower redshifts. 

The impression of luminosity evolution may be mis- 
leading in any case. Spectral spread limits the upper power 
bound of completeness for the RLF in each redshift band. 
At the maximum redshift of the band, radio sources of the 
steepest spectra fall below the survey flux-density limit 
first; the power limit is determined simply from 



268 268 Phni = Sum X X (1 -I- z) 



^S'2.7GHz > Slim, +2.5° >5> -40°, area 3.569 sr. 
tSiiin = 0.25, ^.TGHz > Slim, area 2.278 sr. 

The steps to defining Vohs consist of (1) determining 
P2.7 GHzj the luminosity of the radio source at 2.7 GHz 
(rest frame), and (2) 'moving' the source with its spectrum 
defined by the measured flux densities, from < z < cxd 
to determine in which redshift range(s) it is observable. It 
is observable at a given redshift if (a) its flux density ex- 
ceeds the survey limit 52.7 — 0.25 Jy and (b) its redshifted 
spectrum over the observer's range 2.7 to 5.0 GHz has a 
spectral index > —0.4. We interpolated between measured 
spectral points in the log S^ - log v plane. Despite the rela- 
tively sparse sampling in this plane, combined luminosity 
and spectral effects of 'moving' the source are complex, 
sometimes resulting in a source having two regions of ob- 
servable volume defined by four redshifts. (These effects 
are discussed further in the following section.) In calculat- 
ing the RLF, the contribution of each source is then 

V II iv^ - y'- ) 

/ J I \ * max * mm/ 

i 

where i is usually unity but is sometimes two. Throughout 
the analyses we have used the geometry iJo = 70 km s^^ 
Mpc-i, ntot = 1.0, rJm = 0.3, f^A = 0.7. 

Following these precepts, the radio luminosity func- 
tions calculated for rest-frame powers at 2.7 GHz for 5 
redshift ranges are given in Table |31 



where D is the 'luminosity distance', and amax is the 
minimum (low-frequency) spectral index, i.e. that effec- 
tive index corresponding to the source with 'steepest' ra- 
dio spectrum in the particular redshift range. At lower 
powers within the bin, the RLF will be incomplete for 
such objects, but will remain complete for objects of flat- 
ter spectra. (The limit is well defined for our sample; we 
selected objects of a > —0.4, i.e. the spectral limit was im- 
posed on the 'steep' side, with of course no limit as to how 
'flat' or 'inverted' the spectra might be.) This limit may 
cause RLFs of similar slopes to appear to have a knee at 
similar space densities, mimicking luminosity evolution. 
In previous discussions of space densities it is not clear 
that this limit plus spectral spread have been considered; 
several such studies appear to ascribe a single canonical 
spectrum to every QSO. 

One regrettable result of this power limit is that trac- 
ing the space densities in the higher-redshift ranges down 
to low powers is not possible. Composite RLFs (galax- 
ies plus QSOs) extending over many decades show rel- 
atively few QSOs at low redshifts, where the RLFs are 
dominated b y low-lurninosit y (mostly star-forming) radio 
galaxies ( Sadler et al.ll2002f) . The RLFs of QSOs at high 
redshifts must therefore flatten and drop drastically to- 
wards the lower powers. The dual-population models of 
Jackson & Wall (1999) demonstrate such behaviour. From 
the present data, the limit-lines show that the only con- 
clusion to be drawn is that the RLFs may reduce in slope 
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Table 3. The radio luminosity function p, in units of log(Mpc ^) per Az — 0.5 per A(logP) = 0.4, as derived from 
Sample 2, Table |2| N is the number observed per bin and 1 the mean redshift of the sources in the bin. 
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Fig. 3. The radio luminosity function {Hq = 70 km s"^ Mpc^^, f^tot = 1-0, flm = 0.3, = 0.7) for the QSOs of the 
Parkes 0.25- Jy flat-spectrum sample. The RLF was computed in redshift ranges — 0.5 (red), 0.5—1.0 (green), 1.0 — 2.0 
(blue), 2.0 — 3.0 (orange) and 3.0 — 5.0 (hght blue). Vertical bars show the limits of completeness in power for each 
redshift range, as described in the text. Error bars are 1/VN with N from Table|31 Sky area and redshift-measurement 
completeness have been considered in order to plot true space densities per Az = 0.5, A(logP)= 0.4. 



towards the lower powers. In ^ we show how a differ- 
ent approach can yield some information throughout the 
range of redshifts occupied by the present sample. 

A third presentation of the RLF data is given in Fig.jSJ 
in which space densities are plotted as a function of red- 
shift for 5 ranges of intrinsic power. The initial dramatic 
increase in space density with redshift is evident, with 
densities in the redshift range 1.0 — 2.0 some two orders 
of magnitude above those for objects at redshifts < 0.5. 
Small numbers at the highest redshifts (see Table O and 
the completeness limits at the lower redshifts constrain the 
redshift range observable for each luminosity. In particu- 
lar it is not possible to judge whether the maximum space 
density is a function of radio luminosity. The curves over- 
lap adequately to show self-consistency, and to demon- 
strate the increase in space density from small redshifts 
to z ~ 1.5. Beyond this redshift, the space densities for 



each power range decline, although statistical uncertain- 
ties are substantial. Fig. 13 also indicates such a decline; 
these data therefore suggest a redshift cutoff, at some level 
of significance. 



4. The Redshift Cutoff 

Our preliminary analysis ijShaver et alJll996tl indicated a 
decrease in radio-QSO space density beyond z = 3. Using 
a well-defined sub-sample from the present study, Shaver 
et al. considered the space density of QSOs with P2.7GHZ > 
1.1 X 10^^ W Hz^^ sr^^. On the basis of uniform space 
density, the 25 such radio QSOs seen at z < 4 indicate 
that 15 similar objects would be expected in the range 
5 < z < 7. None was found. From Poisson statistics, the 
difference is significant at the 99.9% level. 
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Fig. 4. The integral radio luminosity function for the 
QSOs of the Parkes 0.25-Jy flat-spectrum sample, com- 
puted in the five redshift ranges of Fig. 01 — 0.5 (red), 
0.5 - 1.0 (green), 1.0 - 2.0 (blue), 2.0 - 3.0 (orange) and 
3.0 — 5.0 (light blue). Vertical lines again indicate limits of 
completeness for each redshift range, due to spectral-index 
spread. The upper curve is the total integrated luminosity 
function, complete for all powers only at the very highest 
luminosities. 
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Fig. 5. Space densities as a function of redshift for 5 power 
ranges, logPa.? = 25.8-26.2, 26.2-26.6, 26.6-27.0, 27.0- 
27.4, and 27.4 — 27.8, denoted by decreasing fine weight. 
Abscissa values are the mean redshifts in each element of 
(AP2.7, Az). For each of the two lowest power ranges, the 
final point with error bar represents incomplete data, as 
these points fall at powers below the spectral cutoffs shown 
in Fig. 121 Upper limits, represented by the single dots, were 
obtained by extrapolating the RLF to this power from 
higher powers. The actual values therefore lie somewhere 
along the two dashed lines. 



This preliminary study drew attention to a possible 
difficulty in the anal ysis due to the curved natu re of some 
of the radio spectra. Ijarvis fc Rawling^ l|2000l) examined 
this in some detail, pointing out the apparently curved na- 
ture of many of the radio spectra involved, and indicating 
how such an effect, a steepening to the high frequencies in 
particular, might reduce or remove the significance of an 
apparent redshift cutoff. Their model-dependent analysis 
used only the highest-power objects and indicated that 
the apparent cutoff on the basis of such objects might 
have a significance level as low as that corresponding to 
2(7. They suggested that establishing the reality of the cut- 
off for such objects to a high level of significance might be 
difficult even with all-sky samples. However, Fig. shows 
that there is no clear majority of sources with spectra 
steepening to the higher frequencies. Moreover, we show 
below that the spectral data in the literature are mislead- 
ing in terms of the proportion of sources showing spectral 
steepening to the higher frequencies. 

Subsequently we have considered alternative methods 
to study space density and redshift distribution, methods 
to utilize the entire sample which can demonstrate sim- 
ple attributes of the space-distribution without recourse 
to modelling the luminosity function or its epoch depen- 
dence. 



4.1. The Power-Volume plane; using the whole sample 

Fig.[3(lcft) shows sources of Sample 3 (TablelSJ in a plot of 
radio luminosity vs. co-moving volume. We need this new 
sample for such a plot. Recall that Sample 1 (Table ^ 
included all sources, not just QSOs, while Sample 2 
(Table EJ, although confined to QSOs, was drawn from 
regions of the survey with different completeness limits. 
In order for a plot of luminosity vs. z (or equivalently, 
co-moving volume) to be interpreted, the sample must 
have a single survey limit. Sample 3 is therefore composed 
of all QSOs from our data table of Paper 1 with survey 
completeness limit at exactly 52.7 GHz = 0.25 Jy (Fig. 1, 
Paper 1), and again at declinations above —40° for reasons 
of radio-spectral completeness. Fig.Qshows lines of survey 
completeness corresponding to 0.25 Jy for three different 
radio spectral indices. 

The plot with co-moving volume on the abscissa rather 
than redshift gives direct indication of space density. 
There is an apparent diminution in the density of points 
at redshifts above ~ 2.5. The question is whether this is 
real and significant. In what follows we test the null hy- 
pothesis that the space density of QSOs at high redshifts 
remains constant and equal to that at 1 < z < 3. 

Redshift information is not complete for Sample 3; in 
order to make comparison with prediction we must esti- 
mate the number of possible objects at 2 > 3. Table O 
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Fig. 7. The Luminosity - Volume plane for the 252 QSOs with measured redshifts, in survey areas with completeness 
hmit <S'2.7GHz — 0.25 Jy (Sample 3, Table|21). Vertical grid lines indicate redshifts, as marked along the top border. The 
curved lines indicate survey completeness limits at 0.25 Jy for three spectral indices. Left: the sources plotted with 
symbols to indicate different spectral indices oi;]^^^^.'- solid dots for —0.4 < a < —0.2; open circles for —0.2 < a < +0.2, 
crosses for -1-0.2 < a < -1-0.6, and stars for the extreme spectral inversions a > +0.6. Rest-frame luminosities (^2.7) 
for this plot were calculated assuming a power-law spectrum given by this index. Note that the dots lie above the 
a = —0.2 limit line; the open circles above the -1-0.2 line; the crosses above the -1-0.6 line; and the stars scatter to 
below the -1-0.6 line. Right: The effects of considering spectral data at frequencies below the Parkes survey frequency 
(2.7 GHz). The values of P2.7 are substantially changed. Crosses indicate the original positions as in the left panel, 
while dots show the positions revised with improved estimates due to incorporation of lower-frequency data. 



presents the summary. The key element is the sub-sample 
of 16 objects in the sample of 268 for which redshifts are 
not available. 

The redshift distribution for QSOs is known to be a 
function of both apparent magnitude and flux density, al- 
beit with huge scatter and only a gentle dependence in 
each case. Thus in order to estimate redshift proportions 
for the objects without such data, we treated the identifi- 
cations made on UKST plates and those from the (deeper) 
CCD observations separately. 

Consider the 11 QSOs without redshifts and identi- 
fied from UKST plates. Of the 242 objects identified on 
UKST plates and with measured redshifts, 8 have z > 3. 
There is no bias in the redshift measurements or lack of; 
and thus we expect 8/242 x 11 = 0.36 of the 11 objects 
to have z > 3. The remaining 5 objects may be treated 
equally; the single unidentified source in the sample (PKS 



0225-065) escaped the CCD identification programme by 
being de-identified later on the basis of an improved radio 
position. Had it been included we can be confident that 
an identification would have been obtained, as it was in 
each of the 87 cases we tried. For these 5 objects, then, 
we use the CCD-identified QSOs with redshifts, totalling 
10 in the sample, for which two redshifts exceeded 3. We 
thus anticipate 2/10 x 5 = 1.0 of the 5 objects will have 
z > 3. The number of objects in the sample with measured 
z > 3 is 10. Thus the number with which to compare pre- 
dictions for z > 3 is 10(observed)-|-1.4(estimated)= 11.4. 
The principal point is that redshift incompleteness does 
not impede our analysis. 

A simple analysis may be carried through on the basis 
of Fig. 13 If we consider QSOs in specific narrow bins of 
luminosity and survey limit imposed by spectral index and 
survey flux limit 5*11111, then such horizontal stripes in Fig. [3 
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Fig. 6. The radio spectra of all sources in Sample 2, 
Table El in their rest frame. Data are at observing fr-e- 
quencies of 0.365, 1.4, 2.7, 5.0, 8.40 and 8.87 GHz, and 
flux densities are normalized by Sq, the interpolated rest- 
frame flux density at 2.7 GHz. The red lines represent the 
brightest sources, those with S'2.7GHz > 2.0 Jy. 

intersecting the curved survey-cutoff lines provide an area 
in the figure in which QSOs can be seen by the survey. 
On the null hypothesis, no redshift diminution, if we now 
split this area into a region with z < 3 and a region with 
2: > 3, we can use the surface density of QSOs in the low- 
redshift area to form an expectation value for the higher 
redshift area. We chose the prediction region to be 1 < 
z < 3 to coincide roughly with the plateau of the 'quasar 
epoch', and we selected the high-redshift region to run out 
to z = 8, the approximate limit to which we could hope 
to see QSOs given our survey limits and the known range 
of luminosity and spectral index. 

This process described above can be refined by reduc- 
ing the stripes of radio power to zero width; each source 
then becomes a predictor, provided of course that the sur- 
vey limit allows it to be seen beyond a redshift of 3. Table^l 
presents results of this analysis under the sub-heading 'sin- 
gle survey cutoff'. The immediate result is the apparent 
one: a prediction of significantly more QSOs at z > 3 than 
the 11.4 'seen'. 

The results reveal a fundamental flaw of this analysis, 
namely what limit line to adopt, corresponding to which 
spectral index. It is apparent from Fig. 13 that adopting 
a = —0.2 is extreme; but even confining the analysis 
to narrow bands of spectral index does not define where 
within that band the survey cutoff or completeness line 
should be placed. The analysis at this point appears to 
confirm what the eye sees in Fig. [3 but shows that taking 
the figure at face value is dangerous. Moreover here we 
have used the 2.7 — 5.0 GHz spectral index, characterizing 
each spectrum as a single power law; spectral curvature 
or indeed any complexity of radio spectrum has not been 
considered. (For low-frequency surveys, the spectral- index 



issue is not so important, because most sources detected 
in them have power-law spectra characterized by an index 
close to —0.75. In corresponding P — z or P — V planes, 
most sources from low-frequency surveys cluster closely 
along or just to the left of the single limit line given by 
this spectral index.) 

The analysis of IShaver et al.l l|l996|) attempted to cir- 
cumvent the difficulties by sticking to powers so high that 
the observational cutoff, the survey completeness limit, 
did not come in to play. In doing so, the available sub- 
sample becomes small and the statistical uncertainties are 
inevitably larger. 

These difficulties suggest the following refinement. 

4.2. Source-by-source analysis: the 'Single-Source 
Survey' 

There is no need to stick to a single survey-limit line in 
the P — V plane. Each source can be considered alone, 
conceptually the result of a survey which found it as a 
single source. For each such 'single-source survey', a limit 
line may be drawn in the plane peculiar to that object and 
incorporating all its radio- spectral information. The pre- 
diction of this object for sources at redshifts above 3 may 
then be added to the predictions from all 'single-source 
surveys' to derive a prediction total. In effect this is us- 
ing the Knax method to predict the number of objects in 
volumes at higher redshift on the hypothesis that space 
density is uniform; it is doing so using the spectral prop- 
erties of each source individually. 

A further advantage in such a process is that there is no 
longer a need to stick to a sample defined by a single flux- 
density limit. To improve statistical weight, all zones of the 
survey can be used, no matter what the flux-density limit, 
provided of course that the value of the 2.7-GHz flux den- 
sity is greater than or equal to the completeness limit for 
the area in which it was detected. (Sources for which this 
is not the case were marked in the data-table of Paper 1.) 
Each source in this analysis contributes a predicted num- 
ber of sources given by the ratio of its accessible co-moving 
volume in the redshift range 3 < z < 8 to that in the range 
1 < z < 3. The sum of all such predicted sources, based 
on all sources observed in the redshift range 1 < z < 3, 
gives us the total number of 3 < z < 8 sources expected 
in the survey for a constant comoving space density. 

A sample appropriate to this analysis is Sample 2 of 
Table 12 giving a total of 379 radio QSOs, 355 with mea- 
sured redshifts. From an analysis analogous to that carried 
out for Sample 3, we estimate that complete identification 
and redshift data would add 1.8 sources to the 16 members 
of this sample observed to have z > 3. 

As a basic analysis of this type, when individual limits 
are applied as described, using the Qf2'7GHz spectral index 
appropriate to each source, a prediction of 51.5 sources in 
the redshift range 3 < z < 8 is obtained (TableQJ, c.f. the 
17.8 sources 'observed'. However a particularly important 
feature of the approach is that it enables incorporation of 



J.V. Wall et al.: The Parkes quarter- Jansky flat-spectrum sample 
Table 4. Predicted numbers of QSOs, 3 < z < 8 



11 



Sample Observed at z 


> 3 Analysis type 


Contributors 


Predictedl" 3 < z < 8 


Single survey cutoff: 








3 11.4 


single cutoff, a = -|-0.2 


73(86) 


64.7(66.7) 


3 11.4 


single cutoff, a = 0.0 


49(57) 


40.3(40.4) 


3 11.4 


single cutoff, a — —0.2 


39(43) 


25.0(25.1) 


Source-by-source analysis, 'Single-source survey': 






2 17.8 


S2.T, S5.0 


71(89) 


51.5(53.5) 


2 17.8 


5o.365, Sl.40, 5*2.7, S5.0 


67(89) 


48.5(54.0) 


2 17.8 


So.365, Sl.40, S2.7, S5.0, S8.40 


53(65) 


28.8(30.1) 


2 17.8 


So.365, Sl.40, S2.7, S5.0, S8.87 Or S8.40 


57(70) 


38.6(38.6) 


2 17.8 


So. 365, Sl.40, S2.7, S5.0, S8.87 


70(88) 


50.5(52.9) 



^calculated in the geometry adopted for this paper; bracketed numbers are for f2m — fltot ~ 1-0, the Einstein - de Sitter Universe 



all radio spectral data. Other investigators (in particular 
iJarvis Sz Rawlingjliooot) have emphasized how important 
the form of the spectrum is. As mentioned in the previous 
section, there are several effects on space density analysis 
and these are illustrated in Figs. [7| and |H1 

1. If the spectrum decreases to low frequencies, perhaps 
having a low-frequency cutoff due to synchrotron or 
free-free absorption, the effective power at the rest- 
frame survey frequency is reduced. The result may be 
a substantially lower position in the P —V plane. This 
reduces the 'headroom' the object has to predict sig- 
nificant contribution at higher redshift; and it may 
remove it from the 'prediction contributors' list en- 
tirely. Significant steepening to the lower frequencies 
of course has the opposite effect, raising its rest-frame 
power, its position in the P — V plane, and increasing 
its prediction. 

We have incorporated spectral data at 1.4 GHz 
from the NVSS survey (Condon ct al. 199^ and at 
0.365 GHz from the Texas Survey ifPouelas et al.ll996(l 
to define the low-frequency spectra for the majority of 
sources in Samples 2 and 3. The results may be seen in 
Fig. 13 The extreme-power objects have convex spec- 
tra and drop down into the pack; but in addition, a 
number of less-luminous objects rise by virtue of steep 
low-frequency spectra. When the prediction is made 
incorporating the low-frequency data (Table 01, these 
effects approximately cancel out and the result differs 
little from the previous estimate: 48.5 objects should 
be seen in the sample at 3 < z < 8, c.f. the 'observed' 
number of 17.8. 

2. A more substantial difference is produced by the in- 
corporation of s pectral data at frequencies higher 
than 5.0 GHz. IShaver et all l(l99fii pointed out 
that high-frequenc y flux dens it ies m easured quasi- 
simultaneously by iGear et all l)l994|) indicated lit- 
tle spectral steepening; and certainly not enough 
to exclude very- high-redshift objects. However the 
iGear et al] l)l994l) sample may not be representative. 



For PKS sour ces a set of flux dens ities at 8.4 GHz was 
measured by IWright et alJ l)l990|) . including a large 
fraction of the sources in both Sample 2 and Sample 3. 
These data suggest that spec tral steepening is more 
common than indicated by the lGear et al.l l)l994(l mea- 
surements, although as Fig.|2|shows, it is not a feature 
of the majority of sources. 

Spectral steepening beyond 5 GHz in the observer 
frame has two effects. First it moves the cutoff line 
upward (Fig. O so that the object in question drops 
from the sample at relatively lower redshifts. Secondly 
when 'moving' the object to some redshift above the 
observed redshift, the spectrum becomes steeper than 
the apparent 2.7 - 5.0-GHz spectral limit of —0.4 
used to define the original sample; 'flat-spectrum' ob- 
jects whose spectra steepen beyond 5.0 GHz (observer 
frame) become undetected as such at higher redshifts. 
To consider the first of these two effects, the cutoff 
(survey limit) line for each source was calculated using 
each 'segment' of the spectrum as redshift is changed. 
This simple interpolation in the log 5*^ - log v plane 
results in segmented cutoff lines for each object in the 
P—V plane as shown in Fig.|S| The second of these two 
effects is illustrated by PKS 2210-257 in the left-most 
panel of Fig.|Hl As the effective spectrum steepens with 
increasing redshift, the accessible region in the redshift 
range 3 < z < 8 is reduced substantially from that 
given by the power cutoff line. The prediction from 
such an object reduces correspondingly. 
There is a third effect due to spectral measurement af- 
fecting the 'observed' region 1 < z < 3, if the spectrum 
steepens at frequencies below the survey frequency. 
The effect is illustrated in the right panel of Fig. |SJ 
objects such as PKS 0017—307 are invisible to us as 
'flat-spectrum' sources at redshifts below a cutoff point 
at which the upward curvature makes them appear 
'steep-spectrum'. This affects only a few sources, but 
for these, it means an increase in prediction contri- 
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Fig. 8. Two sources to illustrate redshift limits in the single-source analysis. For each object, the upper panel shows 
the spectrum in the rest frame, while the lower panel shows the object in the P — V plane. In each case the individual 
cutoff lines are shown as the segmented black curves in the P — V plane, while the smooth coloured curves represent 
completeness limits at 5*2.7 ghz = 0.25 Jy for sources whose spectra are described by single power laws with indices 
as shown. Light grey bars represent the predictive region 1 < ^ < 3, dark grey the region 3 < z < 8 in which the 
object might be visible; the width of these bars is irrelevant. In the case of PKS 2210-257, the individual survey limit 
suggests that the object should be visible out to z=4.54. However, it is not visible to this redshift as an object with 
an observed 2.7 - 5.0 GHz spectrum flatter than —0.4; the spectrum (shown blue) steepens to an effective index of 
< —0.4 at a redshift of 3.22, for which the 'observed' frequencies of 2.7 and 5.0 GHz are shown in the upper panel by 
the red vertical lines. The object therefore drops from the sample at this redshift. PKS 0017-307 does not enter the 
observed region of the diagram as an object with spectral index > —0.4 until a redshift of 1.48 is reached. Again in 
its upper panel, the red vertical hues indicate rest frequencies at 2.7 and 5.0 GHz for the critical redshift of 1.48, at 
which point the spectrum 'flattens' to have a > —0.4, denoted by the blue line. The object runs into its observable 
limit at z = 3.60 in the predictive region as shown. The smaller accessible volume in the 'observed' region results in a 
scaling up of its prediction via the ratio of accessible volumes. 



bution calculated as the ratio of (accessible volume 
3 < 2; < 8) /(accessible volume 1 < z < 3). 

When the available spectral measurements, relatively 
complete for Sample 2 at flve frequencies, are considered 
for each source, the prediction (Table ^ is 28.8 sources, 
differing now from the 'observed' number 17.8 by just 
2.9a. The raised level of the power cutoff does most of the 
damage. It is this use of the data for the Parkes 0.5-Jy sam- 
ple which we believe yields the relativ ely low level of sig- 
nificance for a redshift cutoff found bv lJarvis &: Rawlingj 



However there is a fundamental problem with using the 
8.40-GHz data. This can be shown by using a set of 8.87- 
GHz flux densities measured in 1972 llShimmins fc WaJ 
Il973|) . roughly contemporaneous with the 2.7-GHz sur- 
veys. There are 40 sources in Sample 2 with these 
'old' measurements, one of which, PKS 1532-1-016 at 
Sh s7 h h^ = 1.16 Jy, was not measured bv I Wright et alJ 
lll990|) . There are clearly large flux-density variations at 
8 GHz, the wildest being for PKS 1402-012, 0.67 Jy in 
1972, 0.15 Jy in 1989. If the 1972 8.87-GHz measurements 
are used in preference to the 1989 8.40-GHz measure- 
ments, the prediction is 38.6 QSOs (57 contributors) in 
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the range 3 < 2; < 8. This is very significantly higher 
than the prediction of 28.8 (53) sources using only 8.40- 
GHz data, exceeding the 17.8 sources 'observed' by 4.9cr. 




O 20 40 



predicted number of QSOs at 3 < z < 8 

Fig. 9. The predicted numbers of QSOs at high redshifts. 
The histograms are the results of 2000 trials of end-to- 
end bootstrap testing. The left-most histogram was com- 
piled using high-frequency spectral data at 8.40 GHz and 
8.87 GHz, with the 8.87-GHz (near-contemporary) flux 
densities used in preference if both were available. The 
right histogram represents 2000 results from bootstrap 
testing when only the 8.87-GHz flux densities were used. 
The vertical line to the left indicates the observed number 
of QSOs at high redshifts, as discussed in the text. These 
bootstrap test were carried out in the simple geometry 
J7m — r^tot — 1 but are closely indicative of results for the 
A - dominated geometry adopted here (see Table 0] and 
text). 

The problem is a simple one. Measuring high-frequency 
flux densities some time after the original survey gives 
a biased estimate of the spectrum. Any flux-limited sur- 
vey preferentially selects variable sources in an up-state, 
whereas flux-density measurements many years later re- 
flect sources in a mean state. The result is that the spectra 
are artificially steepened. In the present case the result is 
an underestimate of numbers of objects predicted at high 
redshifts. It is the variations at frequencies above the sur- 
vey frequency which matter in this; variations at the lower 
frequencies are small to insignificant in comparison. 

The result emphasizes how responsive the predictions 
are to flux measurements, and how crucial it is to use con- 
temporaneous measurements. If this much change comes 
about with replacing the 8-GHz flux densities of just 40 
sources with near-contemporary measurements, it is cer- 
tain that the prediction of 38.6 sources based on using 
all the remaining (non-contemporary) 8.40-GHz flux den- 



sities represents an underestimate or lower limit. If the 
8.40-GHz flux densities are ignored and only the 8.87-GHz 
data used as flux densities at frequencies above 5.0 GHz, 
the result is a prediction of 50.5 sources. This must be an 
overestimate. The 8.87-GHz flux densities were measured 
preferentially for bright sources at high frequencies, and 
thus favour objects well above survey-limit lines. We con- 
clude that on the hypothesis of uniform space distribution, 
somewhere between 38.6 and 50.5 sources are predicted to 
have redshifts between 3 and 8 for Sample 2. 

To assess the uncertainties an end-to-end bootstrap ex- 
periment was run for the two possibilities: (i) using as high 
frequency data only the 8.87-GHz (1972) measurements, 
and (ii) using the combination of 8.87 and 8.40-GHz mea- 
surements, with the former taking precedence if measure- 
ments at both fiux densities were available. Because of 
computing time constraints we had to run this experiment 
using the simple geometry of f^tot = = 1- However as 
Table 01 shows, the predictions in this geometry are very 
similar to the predictions of the A - dominated cosmol- 
ogy, the numbers in question being 38.6 and 52.9 for the 
simple geometry vs. 38.6 and 50.5 for the A - dominated 
geometry. The uncertainties should be representative. In 
the bootstrap sampling, random redshifts were assigned 
to each source from the total sample of redshifts. The fiux 
densities for the source were then 'corrected' to that par- 
ticular redshift making use of the measured redshift of the 
object. The results are shown in Fig. [HI Some 2000 trials 
for each of the two possibilities produced no prediction as 
low as the 'observed' number of 17.8 QSOs at z > 3. 

There are two results from this analysis: 

1. The true prediction of numbers of QSOs at 3 < z < 8 
for a uniformly-filled universe lies between 38.6 and 
50.5 objects for Sample 2. This is to be compared to the 
'observed' number of 17.8. An end-to-end bootstrap 
test indicates that the difference is highly significant. 

2. High-frequency flux-density measurements that are 
non-contemporary are dangerous. They bias the spec- 
tral statistics of any variable-flux sample, because sur- 
veys pick out variable sources in their high state, and 
not in their average state. 

5. The form of the evolution 

Figs. 121 and |5l indicate that the form of the evolution, and 
in particular the shape of the decline at high redshifts, 
cannot be inferred directly. As an indirect route, we used 
Sample 2 and proceeded as follows: 

1. Ten redshift limits were set up, from zi — 1.0 to 
zi = 5.5 in steps of Az = 0.5. We then determined the 
combination of P2.7 and effective spectral index yield- 
ing the maximum number of contributors to the RLF 
for the sub-sample complete to each of these redshift 
limits. The numbers of RLF contributors in these com- 
plete sub-samples ranged from 285 at z — 1.0 down to 
just 9 at z = 5.5; the numbers decrease since higher 
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Fig. 10. Above: space density p vs redshift. The individual 
RLFs are each complete from z — 1.0 to z = 5.5 in steps 
of Az = 0.5, in the order purple, brown, orange, dark red, 
light blue, turquoise, blue, green, red, grey. Below: these 
RLFs normalized to agree over the range at 1.0 < z < 2.5. 
The bold black line is a least-squares fit with a polynomial 
of fifth order, given in the text. The grey lines represent 
1000 bootstrap trials. In this process, fits which resulted 
in lines of positive slope beyond z = 5 were rejected. 



luminosities are needed to be complete to the larger 
redshift limits. 

For each of these 10 sub-samples, we computed the 
RLF using the 1/V^nax contributions already calculated 
for previous estimates of luminosity functions. For the 
two samples with redshift limits at 5.0 and 5.5, there 
are no sources in the upper bins, and to use this ob- 
servation we assigned an upper limit of one source to 
each of them. The results are shown in Figure [TUl up- 
per panel. 

Although these RLFs are each now complete from 
z — Q to z — zi, individually they are inadequate to 
trace the whole QSO epoch. Those with smaller values 
of zi by definition cannot reach large redshifts, while 
those with the larger zi are severely noise limited, par- 
ticularly at the low-redshift ends. To combine these 



results to define the overall space behaviour, we nor- 
malized each curve to agree statistically over the range 
1.0 < z < 2.5. We then fitted a least-squares polyno- 
mial through the points, with the results shown in the 
lower panel of Fig. ^| The heavy black line in this 
diagram is 



logp = -12.49+2. 704z-1.145z2 



-0.1796z3-0.01019z^. 



4. Finally a constrained bootstrap experiment of 1000 
runs was used to give an approximation to the un- 
certainty; the grey area in the figure is the result. 

Normalizing as described is only valid if the form of 
evolution is independent of radio luminosity. There is some 
indication in the data (Fig. 1101 upper panel) that the 
turnovers set in at redshifts increasing with luminosity. 
If so, then normalizing as described would broaden the 
maximum, and the overall curve would be representative 
of the overall evolution form, but in no sense formally ac- 
curate. Moreover, the complete individual pieces of lumi- 
nosity function arc not statistically independent, so that 
the shaded area is indicative of uncertainty, but again is 
not accurate in a formal sense. 

These results enable comparison with other high- 
redshift observations. Fig. shows the shaded area of 
Fig. 1101 in the background, with data from recent compi- 
lations of space density as a function of redshift for AGNs 
selected at X-ray and optical wavelengths (left panel) , and 
star-formation rate (SFR) as a function of redshift (right 
panel). 

Agreement with the form of the X-ra y-selecte d QSO 
evolution is remarkably good. ISilverman et all (jl^^) 
found the X-ray decline to agree in form wit h the opti- 
cal decline determined by Fan and co-workers ijFan et alJ 
.2001a c b) from SDSS. Silverma n et al. also showed that 



the COMBO-17 survey resul ts oflWolf et all |2003) follow 
the X-ray data closely. The iHasinger et alJ |20041 X-ray 
AGN results are again in very good agreement with the 
current determination. There is thus general accord be- 
tween the dependence of space density for QSOs found at 
radio, optical and X-ray wavelengths, all showing a rapid 
rise in co-moving space density to z ~ 1.5 followed by de- 
clining space densities at z > 3. However, there are strong 
dependencies of evolution form on luminosi ty, certainly for 
the op tical and X-ray samp l es as noted bv lHasineer et alJ 
l(2004l) and lSilverman et all 1 20041) : and there may be such 
dependence for the current radio-selected sample. The 
dependence on luminosity is well illustrated by the fact 
that the rising curve of (lower-optical-luminosity) QSOs 
selected from the 2dF survey ( Groom ct al. 200^ is dis- 
placed to higher redshifts than the X-ray or radio-selected 
QSOs. The current agreements are illustrative only; analy- 
sis of the significance must await larger samples providing 
better definition of space density evolution as a function 
of luminosity in each wavelength band. 

The relation between QSO space evolution and star- 
formation-rate history is not so clear. The general sim- 
ilarity was first noted in 1997 l|Wal]l llQQsl) : the rise to 
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Fig. 11. Left: Relative space density of QSOs (p) as a function of redshift. The shaded area and black line represent 
the current QSO space-density deteri nination from Figure 1101 Light blue filled circles show the soft X-ray data from 
Chandra and XMM-Newton surveys ("Hasineer et al."2004'), while the dark blue open circles show the results from 
combined Chandra and ROSAT surveys ( Silver man-£t..ai,..20o3) . Space densi t y behaviour of opt icall y-selected QSOs i s 
given by the set of dark red tria ngles: the data are from lSchmidt et alJ ijlQQiT^ . lFan et alJ ( 2001bl) and lFan et alJ 1 20041) . 
The point at z ~ 6 is taken from lFan et al.l l)2004l) and due to conversion between geometries there is uncertainty in the 
ordinate of 0.1. The X-ray and optical QSO data were scaled vertically to match the cur rent determinatio n of space 
density at redshifts 2 to 2.5. Right: Star formation rate density (SFRD; units adopted bv lBlain et al.ll2002() . with the 
shaded area and black line again showing the current estimate of radio QSO space density. Data in optical and ncar^R 
bands are distinguished as squar es: oran ge (Lillv et al. 1996J, g rey l(Co nn olly et al. 1993), green (Steidel et al. 199£ 



light b lue ijBouwens et al.ll2004a|) . black fGiavalisco et alJl2004|) . purple ijBunker et aljr200^ and red l|Bouwens et alJ 
l2004bl) . In general these data have not been c orrected for extinction. Green dots show results of the extinction 
correction of 4.7 suggested bv ISteidel et alJ l)l999 ) for their data: the light blue ba nd represents an estimate of the 
SFR determined from these and other points as analyzed bv lBoiiwens et al l l|2QQ41;i |). Measurem ents from Far-IR and 
sub-mm observations are shown as circles: FIR as orange filled circles from Fl ores et al and sub-mm as dark 

red filled and open circles and blue filled circles from Bouwcns j3t^ aL (2004b), discussed in the text. The current space 
density determination (black curve) was scaled to match the lLillv et al.) |l99^ points, the orange open squares. 



redshifts of ~ 1.5 appears to be of the same form. But 
at this redshift it appears that determinations of the SFR 
from UV, optical and near-IR measures produce a differ- 
ent form of epoch dependence, with an abrupt transition 
to a law almost constant (or diminishing gradually) with 
increasing redshift out to z > 6. There are substantial 
uncertainties in what extinction correction to apply; but 
this form appears to hold whether or not the data are 
extinction-corrected (provided as claimed that the cor- 
rection is not strongly dependent on redshift). The open 
squares of Fig. ^] uncorrected for extinction, show the 
gradual decline, while the band IChapman et alJl2004l) . 



representing a fit to extinction-corrected data, shows a 
star formation rate essentially independent of redshift to 
z > 6. The data from sub-mm o bservations (dark blue 
and red circles, filled and open: IChapman etall l2004|l 
represent estimates from radio-identified sub-mm galaxies 
(red circles), and from these galaxies and sub-mm galax- 
ies combined (blue circles). Chapman et al. point out that 
the similarity of star-formation-rate contributions at z ^ 2 
suggests that the total SFR from all populations may ex- 
ceed the current estimates significantly. 

It remains somewhat puzzling that the sub-mm galaxy 
star-formation-rate appears to drop beyond z = 3, and 
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that it therefore resembles the AGN space-density law 
rather than that of the galaxies detected in the optical 
and near-IR. This may be superficial, in that the points 
are lower limits, and additional components may be found. 
It is perhaps less puzzling that the AGN space-density law 
differs from the overall SFRD in the sense observed. On 
current hierarchical pictures, redshifts beyond 3 represent 
the era of rapid galaxy assembly; and there may be a delay 
before the large enough galaxies have developed to host 
massive black holes, or before the galaxy-building process 
provides orbit organization appropriate to fuel such black 
holes. 

6. Summary 

We su mmarize the results of this and the pr eceding two 
papers l) Jackson et al.ll2002l: iHook et al.ll2003|) . 

1. (Paper 1) The initial goal of the project was to search 
for high-redshift QSOs without dust bias. Optical 
counterparts for essentially all flat-spectrum objects 
in the sample were obtained. No QSOs at redshifts 
greater than 5 were found. 

2. At the fainter flux densities (S'2.7 < 0.4 Jy) and optical 
magnitudes {Bj > 20.0), substantial numbers of flat- 
spectrum radio galaxies are present (Paper 1, Fig. 7). 
These may have influenced previous claims for a hid- 
den popiflation of heavily-reddened QSOs; 'red' QSOs 
ffiieke et al., igYQ.) do not appear to constitute a major 
fraction of the total population sampled here. Amongst 
the reddest of the stellar identiflcations, about one- 
quarter are BL Lac objects, compared to 9 per cent 
with no colour selection, supporting the synchrotron 
interpretatio n for many of the red objects originally 
proposed bv lRieke et al.l |l979) . However, the discov- 
ery of molecul ca.r line emission from a sample of the red- 
der QSOs by iGariin et all l|l99j^ indicates the likely 
presence of dust in some systems, possibly concen- 
trated in dense nuclear tori. 

3. (Paper 2) A composite optical spectrum for flat- 
spectrum QSOs derived from the present sample 
shows clear qualitative differences in comparison with 
radio-quiet composite spectra: the Ly-a(1216A) and 
CIV(1549A) lines are stronger in the current sample. 
There appears to be no significant difference in the 
Ly-a decrement between radio-loud and radio-quiet 
QSOs. 

4. (Present paper) The redshift distribution has been de- 
rived for the quarter- Jansky sample, with considera- 
tion given to the different constituents - red galax- 
ies, BLLac objects an d QSOs. Comparisq i i with the 
space-density models of'Punlo p fc Peacockl ( 19901) and 
Ijackson & Wall (1999) indicates substantial agree- 
ment, although refinement of space-density modelling 
is clearly possible using the current data. 

5. Radio luminosity functions for flat-spectrum QSOs 
have been calculated using the l/V^abs method. These 
show the rapid increase in space density between z = 



and 2 = 1, a flattening between z = 1 and z = 2.5, and 
evidence of declining space density to yet higher red- 
shifts. Attention is drawn to a selection effect produc- 
ing an apparent flattening of the luminosity functions 
towards lower powers in each redshift shell. This effect, 
due to intrinsic spread in radio spectra, may have gone 
unrecognized in previous analyses, leading to overin- 
terpretation of similarities in form of the RLF at each 
redshift. 

6. The reality of a redshift cutoff has been verified using 
the 'single-source survey' Knax method, in which each 
QSO in the sample is used to predict the number of 
sources observable at higher redshifts on the assump- 
tion of uniform space density. The technique is robust 
and model-free, and it makes use of all QSOs in the 
sample, rather than limiting the statistics to the few 
objects of highest radio luminosity. 

7. The result of this analysis has been to demon- 
strate that a space-density diminution exists at z > 
3, at significance levels of > Aa. Precise signifi- 
cance cannot be evaluated, because it is shown that 
spectral data at high frequencies are critical, and 
that non-contemporary flux-density measurements se- 
riously bias the results towards reducing the apparent 
redshift diminution. 

8. In the light of this high-redshift diminution, an esti- 
mate has been made of the overall evolution for the 
radio-loud QSO population. The model-free approach 
has been retained, with the form mapped simply by fit- 
ting a least-squares polynomial to the data. This was 
followed by a bootstrap analysis - 1000 more polyno- 
mial fits - to provide an indication of the uncertainty. 

9. The form of this evolution was compared with that 
determined for X-ray QSOs (Chandra, XMM-Newton 
and ROSAT results) and for optically-selected QSOs 
(primarily SDSS results). Agreement is excellent. It 
was further compared with the evolution of star forma- 
tion rate. Current best estimates of this from optical 
and near-IR data show that although the initial rise 
may be of similar form, there is divergence at redshifts 
beyond 1.5. While the AGN space density dives down 
at z > 3, the SFR appears shows relatively little depen- 
dence on redshift to z > 6. The picture from sub- mm 
measurements is less clear. Uncertainties are large, but 
the SFR determined from sub-mm data appears curi- 
ously to follow the AGN relation more closely than it 
follows the SFR law found from galaxies detected with 
optical and near-IR measurements. 

10. In view of the diminution in space density to high red- 
shifts, radio-loud QSOs would appear to have little role 
to play in the reionization epoch 6 < z < 17. 

The reality of the redshift cutoff for radio-selected 
QSOs and its similarity with that observed for optically- 
selected QSOs leads to the conclusion that the apparent 
cutoff for the latter is real and is not due to dust. This con- 
clusion is_consi^tentwii^ from the CORALS 
survey ((Ellison et al.l200lll2002^ ■ in which a complete sub- 
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sample of radio-selected QSOs from the Parkes quartcr- 
Jansky sample was searched for damped Ly-a (DLA) sys- 
tems. Little significant difference for the comoving mass 
density of neutral gas was found between the estimate 
from the CORALS sample and those from previous (op- 
tically selected) samples. No major population of high- 
column-density absorbers has been missed. Our view of 
the high-redshift diminution of the QSO population does 
not appear to be dominated either by intrinsically dusty 
systems or by a major Universal dust screen. 
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